Electron-beam-induced shift in the apparent position of a pinned vortex in a thin 

superconducting film 
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When an electron beam strikes a superconducting thin film near a pinned vortex, it locally in- 
creases the temperature-dependent London penetration depth and perturbs the circulating super- 
current, thereby distorting the vortex's magnetic field toward the heated spot. This phenomenon 
has been used to visualize vortices pinned in SQUIDs using low-temperature scanning electron 
microscopy. In this paper I develop a quantitative theory to calculate the displacement of the 
vortex-generated magnetic-flux distribution as a function of the distance of the beam spot from the 
vortex core. The results are calculated using four different models for the spatial distribution of the 
thermal power deposited by the electron beam. 



PACS numbers: 74.78.-w, 74.25.Qt, 85.25.Dq 



I. INTRODUCTION 

An important fundamental property of superfluids and 
superconductors is that they admit quantized vortices, 
produced either by rotation of neutral superfluids or by 
applying magnetic fields to superconductors. Numerous 
experimental tools have been used to visualize these 
vortices. In neutral superfluids, charge decoration has 
been used to observe vortices in superfluid helium, 1 and 
ballistic expansion to observe them in Bose-Einstein 
condensates? 2 " In superconducting films, while scanning 
tunneling microscopy^ detects quasiparticles in the 
vortex core, the localized magnetic-field distributions 
generated by the circulating supercurrents open the 
possibility of additional techniques for the observation 
of singly quantized vortices, such as Bitter decoration, 4 



magnetic force 
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microscopy. 
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scanning 



SQUID 

microscopy^ scanning Hall-probe microscopyf&Siifi 
magneto-optical detectionjii Lorentz microscopy^ 
and electron holography^ Low-temperature scanning 
electron microscopy (LTSEM) 14 i 15 i 16 i 17 i 18 i 19 and laser 
scanning microscopj«22*2i have been used to visualize 
vortices in Josephson junctions. 

Recently LTSEM has been used in a new way 22 ' 23 : 24 
to detect the presence of pinned vortices in thin-film su- 
perconducting quantum interference devices (SQUIDs) 
at 77 K. When the scanning electron beam strikes the 
film near a pinned vortex, it locally raises the temper- 
ature, decreases the superfluid density, and increases 
the temperature-dependent London penetration depth 
A(T). The supercurrent circulating around the vortex 
is perturbed, the perturbation being a dipole-like back- 
flow current distribution, which generates a correspond- 
ing magnetic-field perturbation. As a result, the overall 
magnetic-field distribution generated by the vortex is no 
longer centered on the pinned vortex core but is distorted 
in the direction of the beam spot. The shift in the ap- 
parent position of the vortex produces a small change in 
the return magnetic flux threading the SQUID 's central 
hole. As the electron beam rasters across the sample, 
the SQUID output, which is extremely sensitive to such 



flux changes, can be displayed on a video screen. The 
resulting image reveals each pinned vortex as a pair of 
bright and dark spots centered on the vortex, the bright 
(dark) spot corresponding to an increase (decrease) in 
the vortex-generated magnetic flux sensed by the SQUID. 
Since a quantitative description of the above behavior is 
still lacking, I develop in Sec. II a theory to calculate 
the displacement of the vortex-generated magnetic-flux 
distribution as a function of the distance of the beam 
spot from the vortex core. I present results using four 
models for the spatial distribution of the thermal power 
deposited by the electron beam. In Sec. Ill, I summarize 
the theoretical results and discuss possible experiments 
to test the theory. 



II. THEORY 

Consider a vortex pinned at the origin in an infinite 
film of thickness d less than the temperature-dependent 
London penetration depth A(T), such that the relevant 
screening length is the Pearl length 2 ^ A(T) = A 2 jd. Since 
the current density j is nearly constant across the thick- 
ness, we need only consider the sheet-current density 
K = jd. The magnetic induction generated by the vor- 
tex B(p,z) — V x A is described by a vector potential 
A(p, z) that in the plane of the film (z — 0) obeys the 
London equation^ 



A(p,0) 



p AK(p) = -^V 7 



2np 



?Z> (1) 



where K(p) = (2/p, )z x B(p, 0+), O = h/2e is the su- 
perconducting flux quantum, 7 is the phase of the super- 
conducting order parameter, p = xx + yy, p = y 1 x 1 + y 2 , 
p = p/p, and (f> = z x p. When the temperature of the 
film is spatially uniform (T = To), so is the Pearl length 
A = A (T ) = Ao, and the solution of Eq. is well 
known. The sheet-current density isS^Si 



Ko(p) 
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where q = q x x + q y y, q 
z x q. For p <C A , 



Qx + q^ Q= q/q, and <Zi 



solutions are 



K Q (p) w 00o/27r^ o A O /9, 



while for p ^> A , 



Kb(p) ~ #o/tt^oP 2 



(3) 



(4) 



The corresponding magnetic field distribution Bo(p, z) 
is centered on the z axis; thus J B 0z (p 7 z)d 2 p = 4>o 
and / pBo z (p, z)d 2 p = 0. Above the film at a distance 

r = y/ p 2 + z 2 somewhat greater than Ao, -Bo = V x Aq 
appears as if produced by a magnetic monopole at the ori- 
gin; i.e., Bq f» f0o/27rr 2 , where r = p + zz and f = r/r. 

When a sharply focused electron beam scans across the 
superconducting film, depositing thermal power Pq in the 
film and its substrate, this locally raises the film's tem- 
perature around the beam spot at po(t) = xo(t)x+yo(t)y. 
As in Refs. I22U23U24L I consider slow scans such that the 
temperature increment quasistatically follows the elec- 
tron beam and can be written as ST(p — po). This 
temperature increment can be calculated by solving the 
steady-state heat diffusion equation, but this generally 
requires a detailed knowledge of the initial pear-shaped 
spatial distribution of the intensity of thermal energy de- 
position by the e-beam (i.e., the size and shape of the 
interaction volume), the thermal conductivities Kp and 
ks of the superconducting film and substrate, and the co- 
efficient of heat transfer from the film into the substrate 
a s Mil& In the experiments of Ref. l24l done at 77 K, the 
authors estimated that the maximum temperature incre- 
ment AT = ST(0) was only few K. I thus consider only 
small values of AT, such that the deviation of A from 
A , i.e., S A = Ai(p) = (dA/dT)ST(p), is a small pertur- 
bation. What is formally needed in the following theory 



ST(q) 



5T{p)e- lqp d 2 p, 



(5) 



the two-dimensional Fourier transform of the tempera- 
ture increment^ 

To solve for the shift in the apparent position of the 
pinned vortex when AT <C T), it is appropriate to use 
first-order perturbation theory, with A = Ao + Ai and 
A(p, z) — Aq(p, z) + A\(p, z), where quantities with the 
subscript 1 are proportional to AT. Equation Q yields 
the boundary condition 

A 1 (p,0) + paAaK^p) = -pL k x { P - po)K (p). (6) 

It is useful to introduce the Green functions a n (r), 
K{r) = V x a„(r), and k n (p) = (2/p )z x 6 n (p,0+), 
which obey 

a„(p, 0) + p A k n (p) = -5(p)x n , (7) 
where x\ = x, x<i = y, r = (p, z), and n = 1 or 2. The 



/j2 
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1 + 2qA Q 

bw(P! , )= _ /fg 9Mz)+m^) el9 . Pe - 9 \,\ t (9) 



(2tt) 2 1 + 2gA 



kn[P) ~ poJ(2n) 2 l + 2qA 



(10) 



where s(z) = +1 if z > 0, if z = 0, and -1 if z < 0. 
Then Eq. © is solved by 

A^z) = [ D n {p',p )a n (p-p',z)d 2 p\ (11) 

n=l 

2 r 

Bx(p,z) = / D n (p',p )b n (p- p',z)d 2 p\ (12) 

n=l 
2 

K ^P) = E J D n (p',p )k n (p-p')d 2 p', (13) 



where 

2 

D(p,p Q ) = Y D n (p,p )x n = MoAi(p - po)Ko(p). 

n=l 

(14) 

Note that V • fc n (p, z) = and V • K x (p, z) = but that 
V • a„(p, z) 7^ and V ■ Ai(p, z) ^ 0. 

Two important properties of the z component of 
b n (p, z) are that J b nz (p, z)d 2 p = and 



J pb nz (p,z)d 2 p =y^S. 



X m X Z. 



(15) 



Equations and ifTIJl yield / B lz {p, z)d 2 p = and 

J P B lz (p, z)d 2 p = J D(p, p )d 2 p x z. (16) 

Using the above properties of Bo and £?i , we find that 
to first order in AT the electron-beam-induced shift in 
the apparent position of the vortex [i.e., the center of the 
distribution of B z (p, z)] is 



S(Po) 



J pB z (p 1 z)d 2 p _ 1 
fB z (p 7 z)d 2 p 



- [ D(p,p )d 2 pxz. (17) 
■0 J 



Using Eqs. (0), JSJ), and i|14|) and the properties of two- 
dimensional Fourier transforms, we obtain the following 
general expression for the shifts^ 

„, x po dA f°° Sf(q) j . 

s(po) = -3^ / , , o a Mqpo)qdq (18) 



7T dT Jo 1 + 2gA c 

as a function of the be am spot position po = xqx + yoy, 
where po = ^/x\]+y%, po = Po/po, and J n here and in 



3 



later equations is a Bessel function of the first kind of 
order n. Before we can evaluate the integral in Eq. (|18|1 . 
we need to obtain an expression for 8T(q) that provides 
a good description of the experimental conditions. 

Let us consider an experimental situation close to that 
described in Ref. |2J, in which a well-collimated elec- 
tron beam of energy Eg = 10 keV, beam current p 
7 nA, and radius a ~ 5 nm was normally incident upon 
an epitaxial c-axis-oriented superconducting YBa2Cu3 07 
(YBCO) thin film of thickness d = 80 nm on a thick 
SrTiC>3 substrate at 77 K. The range P, i.e., the max- 
imum distance from the point of incidence the electron 
travels until it thermalizes, was estimated to be P w 0.53 
/im. Each incident electron undergoes numerous elastic 
and inelastic scattering processes as it loses energy along 
its diffusion path s, and some of this energy emerges 
from the top surface of the sample as back-scattered elec- 
trons, secondary electrons, x-rays, and electromagnetic 
radiation^ The thermal power delivered to the sample 
is Po = fhEo/e, where the fraction of energy that is 
ultimately converted into heat is of the order of / = 40- 
80% M 

Experiments and Monte Carlo simulations investi- 
gating the electron-scattering and energy-loss processes 
show that the electrons follow relatively straight paths 
when they are incident upon low-Z targets such as car- 
bon or plastic, such that the diffusion cloud resembles a 
paint brushi22iSi On the other hand, when the electrons 
are incident upon high-Z targets such as gold, the more 
frequent large-angle elastic scattering processes cause the 
diffusion cloud to be apple-shapedi^2i2£ The highest den- 
sity of thermal energy deposition is at the center of the 
beam spot, where all the diffusion paths begin. Because 
of the complexity of all the electron scattering and diffu- 
sion processes, it is not possible to describe the density of 
thermal power deposition P v (p, z) with high accuracy by 
a simple analytic formula. 29 However, because a mathe- 
matical expression for P v (p, z) is needed for a calculation 
of 5T(q), which appears in the expression for S(po) in 
Eq. 118|) . I shall present calculations for four different 
approximate models. Each model assumes the electron 
beam is centered on the z axis, P v (p, z) is cylindrically 
symmetric about the z axis, and the volume integral of 
P v (p, z) is equal to Po- 

First, however, it is useful to consider how to solve for 
ST(q) in general for any given density of thermal power 
deposition P v (p,z). At temperatures in the vicinity of 
77 K and above, the temperature perturbation F(p, z) in 
the steady state can be obtained by solving the thermal 
diffusion equation^ — kV 2 P = P v subject to the bound- 
ary condition that dF/dz — at the surface z = 0. This 
assumes that to good approximation the thermal conduc- 
tivities of the superconductor and the substrate are the 
same and equal to k. The general solution is 



where 



F(p,z) = ±- 
Zn 



F(q,z)J Q (qp)qdq, 



(19) 



F(q,z) = — dp'p' I dz'P v (p',z')J (qp') 
K Q Jo 



x (e 



-q\z-z'\ , e -q\z+z' 



I). (20) 



The temperature perturbation at the surface is ST(p) 
F(p : 0) or 



ST(p) 



2n 



6f(q)J (qp)qdq, (21) 



where 



2tt 



oo ,.0 



5T{q) = — I dp'p' dz'P v (p',z')J (qp')e-^'K (22) 
«<7 Jo 



In the limit as q — > 0, ST(q) obeys 



2tt 



8T(q)^— dp'p' dz'P v (p',z') = ^-, (23) 
K <7 Jo J-oo K 1 



and in the limit as p — > oo, 



ST(p) 



Po 

2'KKp 1 



(24) 



where Po is the total thermal power absorbed by the sam- 
ple. For each of the models discussed below, the common 
parameter is P, the maximum range of the electron, and 
it is convenient for later use to express ST(q) and ST(p) in 
terms of the dimensionless auxiliary functions h(u) and 
g(p) via 



Sf(q) = ^h{u), 
Kq 



where u = qR, and 



where p = p/R, 



h(u) = u I J (up)g(p)pdp, 
Jo 

h(0) = 1, 

p OO 

g(p) = / Jo(pu)h(u)du : 
Jo 

poo 

g(0) = / h(u)du, 
Jo 



such that 



ST(0) 



'0 



2itkR 



3(0). 



(25) 

(26) 

(27) 
(28) 
(29) 

(30) 
(31) 



Model KO. To approximate sample heating effects in 
low-Z materials such as carbon, in which the penetrating 
electron beam remains relatively straight, Kanaya and 



4 



Ono^i introduced a model in which thermal power is de- 
posited with uniform density throughout a cylinder of 
radius equal to the the beam radius a and length equal 
to the electron range R, 

Pv,ko(p, z ) = Po/na 2 R, p < a, -R < z < 0, (32) 
= 0, otherwise. (33) 

The function h(u) defined in Eq. H25[) becomes 

Jx{au) (1 - e~ u ) 



(«) 



(cm/2) 



(34) 



where a — a/R, and the steady-state temperature in- 
crease at the surface is given by Eq. (|26|l . where 



9ko{p) 



p'K(k) 



dz' / dp' 

Jo V(P + P') 2 + s ' 2 



(35) 



and K(k) is the complete elliptic integral of the first kind 
of modulus 



^{p + P'Y + ~z 12 

At the center of the beam spot, 
1 



9ko(0) = In 



(VbH?-i) 



(36) 



(37) 



Model C. To approximate sample heating effects in 
higher- Z materials in which the diffusing electron beam 
spreads out radially from the center of the beam spot into 
2tt steradians, one may assume that the thermal power 
is deposited with the following density 

Pv,c(r) = P /27ra 2 i?',r < a, (38) 
= P /2Trr 2 R',a< r < R, (39) 
= 0,r > R, (40) 

where R' = R — 2a/3. This corresponds to the assump- 
tion that the maximum density of deposited thermal 
power occurs at the center of the beam spot but that the 
diffusing electrons lose energy at a constant rate —dE/dr 
as they move radially outward, depositing equal amounts 
of thermal energy into hemispherical shells of increasing 
volume 2irr 2 dr until the electrons thermalize at r = R. 
The function h(u) defined in Eq. 1|25[) becomes 



h c (u) = 1 



1 



-{u[J (u) - (2a 2 /3)J (au)} 



(l-2a/3) 

— [Ji(tt) — (2a/3)Ji(au)) — J2(au)/3u 
-{-ku/2)[J {u)Hx{u) - Jx{u)H Q {u)\ 
+ (nua 2 /3)[Jo(au)Hi(au) — Ji(au)Ho(au)] 
+ [J («) - J (au)]/u}, (41) 

where a = a/R and H n is the Struve function. The 
steady-state temperature increase at the surface is given 



by Eq. where 



9c{P) = 



1 



Mi 



(1 - 2a/3) 

< p < a 

hill 

P 



3a 2 J 



(p-2a/3) 



(1 - 2a/3) 

a < p < 1, 
= 1/P, P>1- 

At the center of the beam spot, 

ln(l/o) + 1/2 



P 



9c(0) 



(l-2a/3) 



(42) 



(43) 
(44) 



(45) 



Model R. Another model, proposed by Reimer, 29 ap- 
proximates the higher density of power dissipation closer 
to the beam spot by assuming that the thermal power 
Pq is deposited uniformly within a hemisphere of radius 
equal to R/2, 



Pv,R{r) 



Po 



[27r(i?/2) 3 /3]' 
= 0,r > R/2. 



<R/2, 



(46) 
(47) 



The function h(u) defined in Eq. 1251) becomes 



h R (u) = 1 - («/2)J (u/2) + Jx(u/2) + (2/u)J 2 («/2) 
+(iru/4)[J (u/2)Hx(u/2) - J 1 (u/2)H (u/2)}. 

(48) 

The steady-state temperature increase at the surface is 
given by Eq. 126(1 . where 



9r(p) 



3-4p 2 , 0<p< 1/2, 
1/P, P > 1/2- 



At the center of the beam spot, 
9r(0) = 3. 



(49) 
(50) 



(51) 



Note that the expression for hn(u) can be obtained from 
that for hc(u) by setting a — 1 on the right-hand side 
of Eq. l(4*T|) and then replacing u = qR by u/2 = qR/2; 
similarly, the result for .9^(0) can be obtained from that 
for gc (0) by setting a — 1 on the right-hand side of 
Eq. H45|) and multiplying by a factor of 2. 

Model B. A similar model, proposed by Bresse^ ap- 
proximates the density of power dissipation by assuming 
that the thermal power Po is deposited uniformly within 
a sphere of radius R/2 centered a distance R/2 below the 
surface, 



Pv,b(p,z) = 



Pn 



[47r(i?/2) 3 /3 
-R < z < 0, 
0, otherwise. 



, vV + {z + R/2) 2 < R/2, 

(52) 
(53) 
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u 

FIG. 1: h(u), the auxiliary function describing the 2D Fourier 
transform 8T(q) via Eq. 1251 . vsu = qR for the four models 
discussed in the text: KO (dotted), C (solid), R(dashed), and 
B(dot-dashed) when a — a/R — 0.01. 
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FIG. 2: g(p), the auxiliary function describing the temper- 
ature perturbation ST(p) via Eq. 12611 . vs p — p/R for the 
four models discussed in the text: KO (dotted), C (solid), 
R(dashed), and B(dot-dashed) when a = a/R = 0.01. 



The function h(u) defined in Eq. I|25[) becomes 



h B (u) 



-u/2 



(54) 



and the steady-state temperature increase at the surface 
is given by Eq. i|26|) . where 



|(P) 



At the center of the beam spot, 
.95 (0) = 2. 



(55) 



(56) 



Figures 1 and 2 show plots of the auxiliary functions 
h(u) vs u = qR and g{p) vs p — p/R for the four models 
KO, C, R, and B in the realistic case that the ratio of the 
e-bcam radius a to the electron range R is a = a/R = 
0.01. 

For each of these models, the electron-beam-induced 
shift in the apparent position of the vortex [see Eqs. i|18|) 
and l|25[) ] can be expressed as 



S(po) = Po 



Pn dA 



G(p /R,2A /R), 



irnR dT 

where G is the dimensionless shift, 

h(u)Jx(p u) 



G(p ,l) 



1 + lu 



du. 



(57) 



(58) 



which is a function of po = po/R and I = 2Aq/R = 
2X 2 (T )/Rd. For models KO and C, G also depends im- 
plicitly upon a = a/R [see Eqs. JjHJl and I0IJ]. 

It is straightforward to numerically evaluate the inte- 
gral in Eq. I|58|l . and Fig. 3 shows plots of G(po,l) vs 



0.15 




FIG. 3: Dimensionless shift G(p ,2) [Eq. (gHJ] vs p = po/R 
for the four models discussed in the text: KO (dotted), 
C(solid), R(dashed), and B(dot-dashed) when a = a/R = 
0.01 and I = 2A /R = 2. As seen in Eq. iJKTJ, this function 
describes the shift in the apparent position of a vortex pinned 
at the origin when the electron beam strikes the film at a 
distance po from the origin. 



Po = Po/R for I = 2Aq/R = 2, which is close to the 
value calculated for the experiments of Ref.|2J. Although 
the curves of G(po,l) vs po are qualitatively similar, the 
distinct differences in shapes and maximum slopes may 
make it possible to determine which of the four models 
best fits the experimental electron-beam-induced shifts 
in vortex position. For po <C 1, the function G(po,l) is 
linear in po; i.e., G(po, I) ~ poG'(l). Since J\{x) ~ x/2 for 
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small x, the initial slope G'(l) can be seen from Eq. 158(1 
to be 



G'(l) 



h(u)u 
1 + lu 



du. 



(59) 



As shown in Fig. 3, the curves for models KO, C, and R 
cross near the origin, and for a = a/ R = 0.01 as p~o — > 0, 
the dotted curve (KO) has the steepest slope G' KO {2) = 
1.25, followed by the solid curve (C) with G' c (2) = 1.07, 
the dashed curve (R) with G' R (2) = 0.53, and the dot- 
dashed curve (B) with G' B {2) = 0.33. As expected from 
the integrands of Eqs. (|58|l and l|59l) . both G and its initial 
slope G' decrease monotonically with increasing values of 
/ = 2Aq/R. For large I, the initial slope approaches the 
value G'(l) = g(0)/2l, and for the four models discussed 
above, this expression gives, when a = 0.01, G' KO (l) — 
2.90/7, G' c (l) = 2.57//, G' R (l) = 1.5/1, and G' B {1) = 1/1. 
The largest values of G and G' occur for very small values 
of I = 2Aq/R. In the limit as I — > 0, the initial slopes are 
G' KO (0) = (l + a-^/l + a 2 )/a 2 or G' KO (0) = 99.5 when 
a = 0.01, G' c (0) = (6-3a)/(12a-8a 2 ) or G" c (0) = 50.1 
when a = 0.01, G' R (0) = 3, and G' B (0) = 2. 

Figure 4 illustrates how the curves of G(po, I) vs po = 
Po/R depend upon I = 2A /R when model KO is used 
with beam radius a = 5 nm, screening length Ao = 0.5 
pm, but electron range R decreasing (with decreasing 
beam energy) from R = 100 pm (I = 0.01) to R = 0.01 
pm (I — 100). Similarly, Fig. 5 (note the different scale 
for the vertical axis) shows how the curves of G(p~o, I) vs 
Pa = Po/R depend upon I = 2Aq/R when model B is 
used with A = 0.5 fim and R decreasing from R = 100 
^m (I = 0.01) to R = 0.01 /xm (I = 100). The general be- 
havior of G(po, I) vs — Po/R is very similar for models 
KO and B, the chief difference being that the slope at 
the origin is much steeper for model KO than for model 
B, as discussed above. Corresponding curves for model 
C would most closely resemble those for model KO, and 
the curves for model R would most closely resemble those 
for model B, as can be seen from Fig. 3 

These results indicate that, for a fixed value of 
Po/ttkR, the largest shifts in the apparent vortex posi- 
tion [Eq. occur for values of I = 2A /R of the order 
of unity or smaller. The order of magnitude of the shift 
when po ~ R is then quickly estimated from Eqs. 10} , 
(Eft, and (3D by noting that Ai - (dA/dT)AT over 
an area of the order of irR 2 and that Kq(R) ~ (f>o/irpoR 2 , 
which yields a maximum shift S(R) ~ (dA/dT)AT, 
where AT ~ Pq/ttkR is the temperature increment at 
the center of the beam spot. On the other hand, for a 
fixed electron-beam current 7b, the maximum shift in the 
apparent vortex position occurs at an intermediate value 
of I, which can be determined for each of the four mod- 
els by noting that the prefactor Pq/ttkR depends upon 
R = 2Aq/1 not only inversely but also implicitly via the 
incident electron energy Eq, as discussed in Sec. III. 
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FIG. 4: Dimensionless shift Gko(po,1) vs po = Po/R, calcu- 
lated using model KO for fixed a and Ao but decreasing R as 
described in the text such that I = 2A /R = 0.01, 0.03, 0.1, 
0.3, 1, 3, 10, 30, and 100. 
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FIG. 5: Dimensionless shift Gb(po,1) vs po = Po/R, calcu- 
lated using model B for fixed Ao but decreasing R such that 
I = 2A /R = 0.01, 0.03, 0.1, 0.3, 1, 3, 10, 30, and 100. 



III. DISCUSSION 

In this paper I considered a vortex, pinned at the ori- 
gin, in a thin film of thickness d when the local heating 
produced by a scanning electron beam focused at po lo- 
cally raises the film's temperature, decreases the super- 
fluid density, and increases the London penetration depth 
A(T). Using first-order perturbation theory, I calculated 
the resulting vortex-generated supercurrent distribution 
and the corresponding distortion of the magnetic-field 
distribution toward the e-beam spot. The resulting ex- 
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pressions for the displacement of the center of the field 
distribution [Eqs. (|57Jl and (|58[l] describe how this shift 
in apparent position depends upon the thermal power 
Po deposited by the e-beam, the range P over which 
this power is delivered, the beam radius a, the ther- 
mal conductivity k, the unperturbed screening length 
(Pearl length) Ao = \ 2 (To)/d, the temperature deriva- 
tive dA/dT of the screening length, and the distance po 
between the e-beam spot and the vortex axis. 

I calculated the shift S(po) using four different models 
(KO, C, R, and B) for the spatial dependence of the 
thermal power deposited by the incident electron beam, 
all models describing the same total thermal power Po 
deposited within P of the center of the beam spot. While 
the results [see Fig. 3] are qualitatively similar for all 
four models, the detailed functional form of S(po) vs po 
is strongly modcl-dcpcndcnt. For example, the slope of 
S(po) vs po at po = is much steeper for models KO 
and C, which account for the greatly increased density of 
power deposition near the center of the beam spot, than 
for models R and B, which assume that the power Po 
is distributed uniformly over much larger volumes. This 
suggests that high-resolution experiments, analyzed with 
the help of the above theory, could be used to determine 
which of the four models gives the best description of the 
thermal power density. 

To observe e-beam shifts in the apparent position of a 
pinned vortex, Doenitz et ali^i used a SQUID with the 
central hole in the shape of a long slot. When a vortex 
was pinned in the body of the SQUID, its return flux was 
measured by the SQUID with high sensitivity, and any 
e-beam- induced shift of the vortex's apparent position to- 
ward or away from the slot resulted in a measurable sig- 
nal, which was displayed as the intensity on a video dis- 
play as the e-beam was rastered across the sample. Such 
an image corresponds to a density plot of S x (xo,yo) vs 
x and y . [Note that S x (-x ,yo) = -S x (x ,yo)-] Shown 
in Fig. 6 is a density plot of S x (xo, yo) calculated using 
model B for the case I = 2Ao/P = 2, roug hly equivalent 
to the experimental conditions of Ref. |2J, where the au- 
thors estimated that R ~ A ~ 0.5 pm. Figure 6 strongly 
resembles the vortex images displayed in Ref. I24L 

Here is an example of a numerical evaluation of 
Eq. H57fl . Assume that the incident electron energy is 
Eo = 10 keV and the beam current is p — 7 nA, as 
in Ref. 124. and that the fraction of the incident elec- 
tron energy that is converted into hea1p2£ is / = 0.6, such 
that P = 42 pW. With a SrTi0 3 thermal conductivity 
k = 18 W/Km at 77 K (Ref. l33l) and an electron range 
R = 0.53 pm (Ref. P /ttkR = 1.4 K, which is the 
temperature perturbation at the center of the beam spot 
according to model B [see Eqs. I|26fl and (|56|l ]. Assum- 
ing A(T) = (140nm)/^i - (T/T c ) 4 and film thickness 
d = 80 nm in A = A 2 (T) /d, as in Ref. one obtains 
A = \ 2 (T a )/d = 0.50 /im and dA/dT = 27.5 nm/K for 
YBCO (T c = 91 K) at T = 77 K, such that I = 2A /P 
= 1.9. For model B, the maximum value of Gb(poJ), 
which occurs at p = Po/R = 1-41, (see Fig. 5) is 0.181, 
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FIG. 6: Density plot of the x component of the vortex 
displacement S x (xo,yo) vs the electron-beam spot position 
(xo I R, yo / R) on the superconducting film for the case that 
I = 2Ao/R = 2, calculated using model B described in the 
text. Positive values of S x are bright, zero values are grey, 
and negative values are dark. The vortex, at the center of 
the plot, shifts to the right (bright areas) when the e-beam 
strikes to the right of the vortex, and to the left (dark areas) 
when the e-beam strikes to the left. The maximum displace- 
ments are shown as the centers of the bright and dark spots. 
A horizontal scan through the center of the plot [S x (xo /R, 0) 
vs xq/R] corresponds to a plot of Gb(po, 2) vs po, as in Fig. 3. 



such that the maximum value of S x (xo,yo) is 7.0 nm at 
xo = 0.75 pm and yo — 0. This corresponds to the center 
of the white spot in Fig. 6. The center of the black spot 
in Fig. 6 corresponds to the value S x — -7.0 nm at xq — 
-0.75 pm and yo = 0. 

To test which of the above four models (KO, C, R, 
or B) is best, similar experiments could be carried out 
to compare S x (xo , 0) vs xo for several values of the e- 
beam radius a for the same incident electron energy and 
accordingly the same range R. According to models KO 
and C, for the case of an electron range R = 0.53 pm, 
as in Ref. 0> there should be a significant reduction in 
the initial slope of S x (xo, 0) vs xo as the beam radius is 
varied from a = 5 nm to 50 nm. According to models R 
and B, however, there should be no significant change in 
slope. 

As another test of the above theory, experiments could 
be carried out over a wide range of electron energies Po 
and corresponding electron ranges R. Assuming that the 
range R oc Pq' 43 , as suggested in Ref.|2!j Doenitz et ali^i 
estimated that P = 0.53 pm at Po = 10 keV. Assum- 
ing that this range-energy relation holds for all energies, 
one finds the results given in Table I when Ao = 0.50 
pm. Comparisons of experimental results with theoret- 
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TABLE I: Calculated values of I = 2ko/R, electron range 
R = 0.53 nm(E /W keV) 143 , and electron-beam energy Eo 
when Aq = 0.50 ^im. 



1 


R (M m ) 


Eo (keV) 


n ni 


i no 


oyu 


0.03 


33 


181 


0.1 


10 


78 


0.3 


3.3 


36 


1 


1 


16 


3 


0.33 


7.2 


10 


0.1 


3.1 


30 


0.033 


1.4 


100 


0.01 


0.62 



ical plots of S x (xq,0) vs Xq obtained from calculations 
such as those shown in Figs. 4 and 5 would provide a 
stringent test of the four models. 

Experiments done at various temperatures above 77 K 
would provide a further test of the theory. In Ref. Izil 
Doenitz et al. estimated the value of Ao = \ 2 (To)/d w 
0.50 fim for YBCO (T c = 91 K) of thickness 80 nm at 
To = 77 K by assuming A(T) = (140 ma)/y/l- (T/T c ) 4 . 
With R = 0.53 fim, this gives a value oil = 2Aq/R = 1.9. 
Increasing the temperature would not only increase the 
value of I and decrease the magnitude of G (see Figs. 4 
and 5) but also increase the magnitude of dk/dT. For 
large I, the combination of these two competing effects 
would lead to an overall increase in the magnitude of 
S x (xo,yo) oc (l/Ao)(dA/dT) according to Eqs. and 
(|58|l . With the temperature T approaching T c , however, 
it is likely that the vortex would become depinned and 
follow the warmer beam spot, such that the above theory, 
which assumes that the vortex remains pinned as the 
electron beam scans over it, would no longer apply. 

Reducing the temperature to lower values would de- 
crease the value of I, thereby increasing the magnitude of 
G. On the other hand, it is likely that the reduced value 
of dk/dT would more than compensate for this effect and 



lead to an overall decrease in the magnitude of S x (xo, yo). 
Going to low temperatures, however, would put the ex- 
periments in a temperature range where the above theory 
for the temperature increment involving only the thermal 
conductivity K is no longer valid. 18 

It is important to note that the above theory is ex- 
pected to be valid for experiments on superconduct- 
ing films at liquid-nitrogen temperatures or above. At 
much lower temperatures, including liquid-helium tem- 
peratures, the above calculations of the temperature per- 
turbation [Eqs. ©-ED] are no longer valid, because it 
is then essential to take into account the thermal bound- 
ary resistance between the superconducting film and the 
substrateii^i^ii 8 ^ In such a case, the theory involves an 
additional length scale, r\ — (nd/ds) 1 / 2 , the thermal heal- 
ing length, where a s is the coefficient of heat transfer 
from the superconducting film into the substrate 

Although the results given in this paper apply strictly 
only to a vortex in an infinite thin film, I expect that 
when R > A , Eqs. lfST|) and i(5*S)l also apply to good 
approximation to a vortex in a thin film with finite lateral 
dimensions, provided that the distance of the vortex from 
the edge of the film is more than a few R. This is because 
the current-density perturbation is largest within an area 
of order R 2 when the electron beam is roughly a distance 
R from the vortex. However, when the distance of the 
vortex from the edge of the film is approximately R or 
smaller, the above calculation would have to be redone, 
taking into account how the supercurrent is modified near 
the film edge. 
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